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An explanation for the curious mass loss history of massive stars: 
from OB stars, through Luminous Blue Variables to Wolf-Rayet 

stars 
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Abstract. The stellar winds of massive stars show large changes in mass-loss rates and terminal velocities during 
their evolution from O-star through the Luminous Blue Variable phase to the Wolf-Rayet phase. The luminosity 
remains approximately unchanged during these phases. These large changes in wind properties are explained in 
the context of the radiation driven wind theory, of which we consider four different models. They are due to 
the evolutionary changes in radius, gravity and surface composition and to the change from optically thin (in 
continuum) line driven winds to optically thick radiation driven winds. 
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1. Introduction 

In this paper we present an explanation for the drastic 
changes in mass loss rate and wind velocity during the 
evolution of massive stars. 

The stellar winds of massive stars show a large vari- 
ation in mass-loss rates and terminal velocities during 
their evolution from O-star through the Luminous Blue 
Variable (LBV) phase to the Wolf-Rayet (WR) phase. O- 
stars have a relatively small mass-loss rate (M ~ 10 -6 
to 10~ 5 Mq yr _1 ) but high wind velocity (vqo ~ few 
10 3 km s _1 ), LBVs have a high mass-loss rate (few 10 -5 
Mq yr _1 ) but a small wind velocity (few 10 2 km s _1 ) and 
WR-stars have a high mass-loss rate (few 10 -5 M© yr -1 ) 
and a high wind velocity (few 10 3 km s _1 ). The luminos- 
ity remains approximately unchanged during these phases. 
The question is: what causes these changes in mass-loss 
rate and in wind velocity? 

The mass-loss rates of OB-stars and LBVs during their 
quiescent phase have been explained in terms of radiation 
driven winds, where the driving is done by multitudes of 
spectral lin es (C astor et al. 1975 , Pauldrach et al. 1986 , 
Vink et al. E000D . Line driven winds are optically thin in 



their continuum. Nugis & Lamers ( 2002 , hereafter NL) 
have shown that the mass-loss rates of Wolf-Rayet stars 
can in principle be explained by optically thick radiation 
driven winds. In this paper we explain the large changes 



in M and in terms of transitions between four types 
of radiation driven wind models, due to changes mainly 
in the stellar parameters and to a smaller extent also in 
the surface abundance during the stellar evolution. We 
discuss the properties of two types of optically thin and 
two types of optically thick radiation driven wind models, 
which we then apply to nine characteristic massive stars 
with increasing evolution stages. 

2. Evolutionary changes of the stellar parameters 
of massive stars 

The evolution of stars initially more massive than abou t 
50 Mq proceeds as follows (e.g. Maeder & Meynet 1987 ). 
After the main sequence phase the star expands and be- 
comes a blue supergiant with a radius of about 10 2 Rq and 
with an enhanced He and N surface abundance. For a rea- 
son that is poorly understood the star becomes an unsta- 
ble LBV with multiple outbursts (Humphreys & Davidson 
1994; Leitherer 1997). Even during quiesence the mass- 
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loss rates of LBVs are significantly higher than during the 
main sequence phase but the wind velocity is much lower. 
After the LBV-phase the star contracts and becomes a 
WR star with a N-rich (WN-type) and later a C,0-rich 
(WC-type) surface. The hydrostatic radius of the star is 
only a few Rq, but the high mass- loss rate produces an 
optically thick wind with a sonic point radii of about 15 
to 30 Rq for hydrogen-rich WNL-phase and about 1.5 to 
10 Rq for the H-poor WNE/WC phase (see NL). Very 
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massive stars (M > 60 Mq) may skip the LBV phase 
(Bohannan & Crowther 1999). 

We adopt nine typical massive stars of about the same 
luminosity and initial mass to characterize the changes 
in the stellar and wind parameters during the evolution 
of a massive star. They represent different evolutionary 
phases. The stars are listed in Table [t] in order of increas- 
ing evolutionary stage: one Of-star (£ Pup), one LBV (P 
Cyg) , three stars with spectra in between those of Of and 
late-WN (HD151804, HD 152408, HD152386), three WN- 
stars (WR 105, WR 136 and WR 139) and one WC-star 
(WR 111). Notice the large changes in M and and 
the general increase of the atmospheric He/H-ratio and of 
the momentum transfer efficiency, rj = Mvoo/{L/c), with 
evolution stage. The radii of the four genuine WN and 
WC-stars in this table are the radii of the sonic point, 
derived by NL, and the value of T e g is at that radius. 

3. Radiation driven wind models 

3.1. Optically thin line driven winds 

For line driven wind models, which are optically thin in 
the continuum, the predicted terminal velocity is 



Voc ^ Cf d \/a/(l - a) 



(1) 



(Castor et al. 1975; Kudritzki et al. 1989] ), where a is a 
force multiplier parameter with a ~ 0.5 to 0.7 for hot 
massive stars of T c ff <^ 8000 K and v csc is the effec- 
tive escape velocity, i.e. corrected for radiation pressure 



by electron scattering. Lamers et al. (1995) have shown 
that C fdv /a/(l -a) ~ 2.7 if T off > 21 000 K and 1.3 if 
10 000 < T cff < 21 000 K. Observations and theory both 
show that the mass-loss rate of a line driven wind increases 
by about a factor 5 and the terminal velocity decreases by 
about a factor two when the effective temperature of the 



star drops below about 21 000 K (Lamers et al. 1995 ; Vink 
et al. 1999). This is the bi-stability jump, which is due to 
the change in ionization in the lower wind layers near the 



sonic point (Vink et al. 1999). Detailed calculations of line 



driven wind models with multiple scattering by Vink et al. 
( 2001 ) have shown that the mass-loss rate of galactic OB- 



stars in Mq yr is 

logM = -6.86 + 2.194 log(L/10 5 ) - 1.313 log(M/30) 
+0.933 log(T eff /40 000) - 10.92{log(T cff /40 000)} 2 (2) 

on the "hot side" of the jump, T cff > 21 000 K, and 



logM = -6.39 + 2.210 log(L/10 5 ) 
+ 1.07 log(T cff /20 000) 



1.339 log(A//30) 



(3) 



on the "cool side" of the jump, 10 000 < T cff < 21 000 K, 
with M and L in solar units. We will use Eqs. |l], || and [| 
to predict Voa and M for line driven winds. 



3.2. Optically thick radiation driven winds 

NL have shown that in optically thick radiation driven 
winds the opacity has to increase outwards at the sonic 
point. They showed that 
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M 



4.66 x 10 



-29 



M 



'(1+7) 



(4) 



in Mq yr -1 , where R s (in Rq) and T s (in K) are the radius 
and temperature at the sonic point, M is in Mq, v s is the 
isothermal sound speed, 7 is the mean number of free elec- 
trons per atom and \x is the mean atomic weight in atomic 
mass units. We see that M of an optically thick radia- 
tion driven wind is proportional to T s . At first sight this 
might suggest that an arbitrary high mass-loss rate can 
be reached by moving the sonic point deeper into the star 
where the optical depth and the temperature are higher. 
However, this is not the case, because the transition from 
subsonic to supersonic velocity at the sonic point sets re- 
quirements for the opacity and its gradient. The opacity 
and its gradient at the sonic point are (see NL) 



A-kcGM 



dx 
dr 
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(5) 
(6) 



where a\ = k("f + l)/(/im u ). Eqs. || and |^ imply that the 
transonic transition can only occur in the layers where the 
opacity increases outwards and where it reaches a value 
set by the luminosity and mass of the star ( Eq. |5 ). ^,From 
the OPAL-opacity tables (Iglesias & Rogers [1996) we find 
that this occurs only in limited temperature regimes where 
x(T) shows a bump. These regimes are in the ranges of 
156 000 <T S < 162 000 K, and 37 000 < T s < 71000 K, 
where respectively a large and a small Fe-opacity peak oc- 
cur (see NL) . We will use Eq. || to predict M for optically 
thick winds. For we adopt the scaling predicted by the 
models. We derived from WR models of NL that 



(2 + 0.5K 



(7) 



with v% sc at the sonic point. We will use this scaling law 
to estimate Voo for optically thick radiation driven winds. 
Notice that increases with decreasing sonic radius. 

3.3. Four types of wind models 

The description above has shown that radiation driven 
winds from hot stars come in four types: 

(1) line driven winds which are optically thin in the con- 
tinuum for stars with T c g > 21 000 K, i.e. on the hot 
side of the bi-stability jump: "line hot" models. For these 
winds we adopt Eqs. ^ and ||. 

(2) line driven winds for stars with T e g < 21 000 K, i.e. on 
the cool side of the bi-stability jump: "line cool" models. 
We adopt the Eqs. [I] and |^ for these models. 
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Name 


Type 


log L 


R w 


T e ff 


M (2) 


v 

Vesc 


N He /N H 


M 






V 


Ref (4) 






L 


Rq 


K 


M Q 


km s~ 




M Q yr 


— i 


km s~ 






C Pup 


04 If 


6.00 


19 


42 000 


70 


953 


0.15 


3 10" 


(i 


2250 


0.32 


PU 


P Cyg 


B1.5 Ia+ 


5.86 


76 


19 300 


23 


223 


0.30 


2 10 _ 


5 


210 


0.28 


PP 


HD151804 


08 Iaf 


5.84 


37 


26 700 


46 


581 


0.25 


1.2 10 


-5 


1445 


1.24 


CB 


HD152408 


WN9ha 


5.80 


32 


27 600 


44 


648 


0.67 


2.4 10 


-5 


995 


1.79 


CB 


HD152386 


WN9ha 


5.82 


33 


27 000 


46 


624 


0.27 


2.7 10 


-5 


1650 


3.45 


BC 


WR 105 


WN 9 


5.81 


26 


32 100 


22 


412 


0.44 


2.8 10' 


-5 


1200 


2.6 


NL 


WR 136 


WN 6b 


5.73 


4.6 


73 000 


19 


900 


1.9 


6.3 10 


-5 


1600 


9.2 


NL 


WR 139 


WN 5 


5.21 


2.0 


82 000 


9.3 


1129 


5.0 


0.9 10' 


-5 


1800 


4.9 


NL 


WR 111 


WC 5 


5.31 


2.3 


81 000 


10.6 


1160 


oo 


1.0 10' 


-5 


2415 


5.8 


NL 



^ For WR-stars the mean sonic radius of the models Al and Bl of NL with the corresponding T fs = (L/ (Ana R 2 )) ' 25 is listed. 
Masses of O and WN9ha-stars were derived from evolutionary models. 



(3) 



is the effective escape velocity at the radius R, corrected for radiation pressure by electron scattering. 



( 4) PU = Puis et al. (1996); PP = Pa uldrach & Puis (1990); CB = Crowther & Bohannan (1997); BC = Bohannan & Crowther 



(1999); NL = Nugis & Lamers (2002) 



(3) "thick cool" continuum driven winds, with the sonic 
point in the temperature range of 38 000 < T s < 71 000 
K, where the small opacity bump occurs (see NL). For 
these stars we adopt Eqs. [| and with T s = 40 000 and 
70 000 K. 

(4) "thick hot" continuum driven winds. For these stars 
we adopt Eqs. | and @ with T s ~ 160 000 K. 

4. Predicted radiation driven mass-loss rates and 
velocities 

We apply the predictions of the radiation driven wind 
models to the nine stars. The resulting values of M and 
Voc are listed in Table 2 for four models: "line cool" or 
"line hot" , "thick cool" with T s = 40 000 and 70 000 K, 
and "thick hot" with T s = 160 000 K. For the line driven 
wind model of P Cyg we adopt the predictions for the cool 
side of the bi-stability jump ("line cool") whereas for the 
other stars we adopt the "line hot" models. 

Comparing the predicted values of M and of the 
four models with the observed values, we can determine 
which model fits best. In this comparison we have given 
more weight to the mass-loss rate than to the terminal ve- 
locity, because the predicted Voo of the line driven models 
depends on the mass of the star, which is not well known, 
and of the optically thick winds is not well predicted 
by NL models. For the star £ Pup the models "line hot" 
and "thick cool" with T s = 40 kK predict almost the same 
mass-loss rates. We adopt the "line hot" model, because 
this star is considered to be the prototype of a line driven 
wind (e.g. Pauldrach et al. 1994). In the case of the 08 Iaf 
star HD 151804 M is in between the predicted values of 
the "line hot" and "thick cool" (40 kK) models. The ob- 
served value of u oo agrees better with the "line hot" model. 
The wind of this star may be of intermediate type. For the 
other stars the choice of the best fitting model is quite ob- 
vious (last column of Table ||). Notice that for the "thick 
cool" models the observed values of M agree better with 
those predicted for T s ~ 40 kK than for T s ~ 70 kK. 



5. Discussion and conclusions 

^From the comparison between the predicted and the ob- 
served values of M and Vqo , we can explain the changes in 
mass-loss rate and wind velocities during the evolution of 
massive stars in terms of the four wind models. 

1. O-stars on the main sequence and shortly thereafter 
have winds driven by lines on the hot side of the bi- 
stability jump. 

2. The winds of LBVs are radiation driven by lines. The 
increase in M and the decrease in Woo from O-star to 
LBV is due to the crossing of the bi-stability jump. 
(However, not all LBVs become cool enough to reach 
the bi-stability jump: Leitherer 1997; Lamers 1997) 

3. When the star has lost sufficient mass and the atmo- 
sphere has been He-enriched sufficiently to contract to 
the WNL-phase, the wind becomes optically thick and 
the sonic point moves into the region where the contin- 
uum opacity shows a small bump so that it can initiate 
an optically thick wind. This results in a (small) in- 
crease of M and a large increase in i>oo ■ The character- 
istics of the star HD 151804 suggest that the transition 
from line driven winds to "thick cool" is gradual. 

4. When the star looses more mass and the surface be- 
comes H-poor its wind may either stay "thick cool" (at 
T s ~ 40 kK) or become "thick hot". In the first case 
the star may appear as a H-poor WN7 or WN8-star 
(not studied here). In the latter case the star appears 
as a WNE-star, similar to WR 136 and WR 139. In 
both cases the star has a high M and a high i>oo- The 
transition from a "thick cool" to "thick hot" wind is 
determined by the variation of Xs (see Eq. ^|). During 
the WNE-phase the M/L-ratio increases (Schaerer & 
Maeder 1992). When Xs increases to values above ~0.5 
cm 2 g _1 , the sonic point has to move to high tempera- 
ture regime, because such a high value of Xs is reached 
only near the main iron opacity peak (NL). When the 
star evolves directly from Of to WNLh both the M/L- 
ratio and Xs decrease. 
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Table 2. Predictions for line driven (optically thin continuum) and optically thick winds versus observations. 



Star 


Tvne 






log M 


(Mn vr'l 

V JW j 1 7 






(km s _1 ) 




Best 






thin 


thick 


thick 


thick 


otas 


thin 


thick 


obs 


in nil pi 

111UUC1 






line 


cool 
40 kK 


cool 
70 kK 


hot 
160 kK 




line 








(, rup 


04 If 


-5.12 


^ ^4 


-4.44 


9 


^9 






2250 


lino ri i~il~ 
1111c 11L1L 


P Cyg 


Bl Ia + 


-4.35 


-3.28 


-2.19 


-0.58 


-4.70 


290 


446 


210 


line cool 


nn i p: i snzi 

ilL' _LtJ-LOU4: 


OS Tnf 


t\ 71 


A M 


-0.4:^ 


1 si 


A Q9 


ijuy 


1 1 R9 


1 AA^ 


line / LiiieK. ; 


HD 152408 


WN9ha 


-5.68 


-4.73 


-3.65 


-2.02 


-4.62 


1750 


1296 


955 


thick cool 


HD 152386 


WN9ha 


-5.78 


-4.67 


-3.57 


-1.96 


-4.56 


1685 


1248 


1650 


thick cool 


WR 105 


WN 9 


-5.09 


-4.68 


-3.59 


-1.97 


-4.55 


1037 


824 


1200 


thick cool 


WR 136 


WN 6b 


-5.50 


-6.94 


-5.84 


-4.22 


-4.20 


2430 


1800 


1600 


thick hot 


WR 139 


WN 5 


-6.50 


-7.73 


-6.64 


-5.03 


-5.05 


3048 


2258 


1800 


thick hot 


WR 111 


WC 5 


-6.33 


-7.70 


-6.60 


-4.99 


-5.00 


3132 


2320 


2415 


thick hot 



5. It is difficult to predict the dependence of M on L for 
optically thick wind models accurately because R s and 
hence also T e g(R s ) is not well known. The wind mod- 
els of NL for WNE/WCE stars predict that R s oc L 7 
and because M oc L 6 (Schaerer & Maeder |l992|) , it 
follows that M oc R s 3 /M oc L 1 - 5 (Eq. |). This depen- 
dence agrees well wi th th e empirical relation derived 
by Nugis & Lamers fl2000|) for WN stars. For the WC 
stars Nugis & Lamers (2000) found empirically that 
M oc _L 84 but with strong dependence on chemical 
composition. On the other hand, in a recent study of 
LMC WC-stars Crowther et al. (A&A in press) found 
a strong dependence of M on L (M oc L 1 - 38 ) which 
agrees well with our predicted dependence for optically 
"thick hot" wind models. 

6. We find that the changes in M and Vco during the 
evolution of massive stars are mainly due to changes in 
the stellar parameters and to a lesser degree to changes 
in the surface composition. 

We have shown that the changes in M and Voo during 
the evolution of the massive stars from O-star to WN-star 
are due to the adjustment of the wind to the changing con- 
ditions, mainly the M , L, R and surface composition. We 
have not explained "how" these changes occur. The tran- 
sition from line driven winds to "thick cool" winds is most 
likely due to the formation of a bump in the opacity curve 
when the He/H ratio increases. This transition can occur 
gradually because the temperature of the sonic point of a 
"line hot" wind overlaps with the range for "thick cool" 
winds. The transition from "thick cool" to "thick hot" is 
probably due to the fact that the hydrostatic radius of a 
WR-star shrinks and the M/L-ratio increases when the lu - 
minosity of the star decreases (Schaerer & Maeder 1992| ). 
This transition cannot be gradual, because the sonic point 
temperature in these models is very different, so the wind 
must be restructured during this transition. 

We point out that at present, the optically thick wind 
models do not provide accurate predictions for v^, so 
we adopted an empirical scaling law. However, Schmutz 



supersonic part of the winds of WR-stars to explain the 
observed high values of v^. 
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